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ABSTRACT 

Using a one-dimensional height integrated model, we calculate the evolution of an 
unequal mass binary black hole with a coplanar gas disc that contains a gap due to the 
presence of the secondary black hole. Viscous evolution of the outer circumbinary disc 
initially hardens the binary, while the inner disc drains onto the primary (central) black 
hole. As long as the inner disc remains cool and thin at low Mcxt (rather than becoming 
hot and geometrically thick) , the mass of the inner disc reaches an asymptotic mass 
typically ~ 10"'^ — 10~^Mq. Once the semimajor axis shrinks below a critical value, 
angular momentum losses from gravitational waves dominate over viscous transport 
in hardening the binary. The inner disc then no longer responds viscously to the 
inspiraling black holes. Instead, tidal interactions with the secondary rapidly drive 
the inner disc into the primary. Tidal and viscous dissipation in the inner disc lead 
to a late time brightening in luminosity, L oc ij^jjj^g, where tminus is the time prior to 
the final merger. This late time brightening peaks ~ 1 day prior to the final merger 
at ~ O.lLEdd- This behavior is relatively robust because of self regulation in the 
coupled viscous-gravitational evolution of such binary systems. It constitutes a unique 
electromagnetic signature of a binary supermassive black hole merger and may allow 
the host galaxy to be identified if used in conjunction with the Laser Interferometric 
Space Antenna (LISA) localization. 

Key words: black hole physics ~ accretion, accretion discs - binaries: general - 
gravitational waves - galaxies: active - galaxies: nuclei ~ quasars: general 



1 INTRODUCTION 

There is now strong evidence that supermassive black 
holes (SMBHs) inhabit the spheroids of most galaxies 
l|Magorrian et al.lll998l : lTremaine et a l. 2002). Galaxy merg- 
ers, which are a natural consequence of hierarchical assembly 
of galaxies, should cause the SMBHs in each galaxy to form a 
binar y in the newly assembled host galaxy ( Begelman et al.l 
Il980h . These binary black holes (BBHs) are believed to coa- 
lesce due to an uncertain combination of dynamical friction, 
three body interactions, interaction with viscous circumbi- 
nary discs, and gravitational wave (GW) losses. Observa- 
tions of the GW that these systems produce is one of the 
primary scientific objectives of the proposed Laser Interfer- 
ometer Space Antenna (LISA). In addition to their useful- 



ness as probes of general relativity in the strong field regime, 
GWs can also be used to measure the luminosity distance 
to hig h accuracy, i.e., SDh/Dh at the per cent level in many 
cases (|Hughesll2003 : lHolz fc Hughes^l2005^ . 

The high accuracy that these "standard sirens" of- 
fer would allow these GW sources to constrain cosmologi- 
cal models if the host galaxy can be identified, so that a 
source redshift can be measured. By itself, LISA observa- 
tions will provide sky localization er rors that are typical! 



few X 10 deg (see for instance Holz fc Hughei 



i^ pically 
il2005l: 



i Kocsis et all |2007| . l2008l : iLang fc Hughej|2008l : lArun et all 
|200^ . in which there may be ~ 10* galaxies Q Hence, much 
recent work has focused on finding an una mbiguous elec - 
tromagnetic counterpart to a BBH merger ISchutzl \l9M ) 
to help identify its host galaxy. In addition, the known BH 
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^ The 3-D localiz ation of LISA would g i ve ~ 10 candidate galax- 
ies per arcmin-^ i IHoIz fc Huehesll2005i ; iKocsis et al.ll2006l) . 
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masses and spins from GWs, an observed electromagnetic 
counterpart would also provide strong constraints on mod- 
els of accretion discs and quasars/ AGNs. 

At present, most studies of these electromagnetic 
counterparts to BBH mergers have focused on the post- 
merger phase. At the time of merger, a fraction of the 
rest mass of the two BHs is radiated as GW and the 
merged BH recei v es a kick between 100 — 4000 km s~^ 
(Bckonstein '1973; Favata ct al j l2004l: 
Herrmann et al. 2007; Go nzalez et alJ 
20081 '). If a circumbiiiarv gas disc surrounded the two BHs 
(pre-merger), the sudden decrease in BH mass and the kick 
provided by the BH may produce an observable electro- 
magn eti c signature jBode fc PhinnevI 12007 : Lippai et al. 
20081 : ISchnittman fc Krohkl l2008l : IShields fc Bonnine 



Baker et al.) l2006l: 



20071 : iBovle et al I 



20081 : lO'Neill et al.l l2008l : iBode fc PhinnevI l2009l l. On a 
longer timescale, the circumbinary disc would accrete 
into the newly merged BH, leading to AGN activity 



Milosavlievic fc P hinncv 2005') or perhaps an offset quasar 



Madau fc Quataert, 2004 ). 



There are potential difficulties associated with iden- 
tifying these electromagnetic signatures of BBH merg- 
ers. The signatures associated with kicks are likely to 
be weak and difficult to distinguish from AGN variabil- 
ity, if th ey are a few percent of the background disc 
emission (ISchnittman fc KrolikI [ioosi : iKocsis fc Loebll2008l : 
lO'Neill et al.ll2008D . Delaved AGN and offset AGN activity, 
while significantly more prom ising, might have a substantial 
delay l|Milosavlievic fc Phinnev 2005 ; Loeb 2007; ) . 

The subject of this work is to explore elec- 
tromagnetic emission ass ociated with the pre-merger 
phase. Various autho r s llMilosavlievic fc Phinnev 20051: 



Armitage fc NataraianI l2005l : iMacFadven fc Milosavlievid 
2008^) have studied the physics of such viscously driven 
BBH mergers. These studies suggest that the interaction 
between the disc and the BBH s can induce a small eccentric- 
ity in the disc and the binary llArmitage fc Nataraianll2005l : 
IMacFadven fc Milosavlievid l2008l : ICuadra et al.1 120091 ). In 
addition, some of this gas may form accretion discs around 
the t wo BHs leading to binary AGN activity ijPotti et al.l 
I2OO7I ). 

In this paper, we build on the work of 
lArmitage fc NataraianI l|2002l ) and show that although 
the inner disc drains when the binary is well separated, a 
small amount of remnant or fossil gas can remain around 
the primary BH, giving rise to a distinct electromag- 
netic signature. The amount of fossil gas is very small 
(~ 10-^Mo), but even this tiny amount of gas is sufficient 



to drive the luminosity of the binary to 
before merger. 



0.1 1/Edd one day 



The plan of this work is as follows. In we summarize 
the basic scenario. We then present the relevant physics in 
^ In 21 we argue that a small amount of fossil gas could 
exist in the inner disc between the primary and secondary 
BHs and that this fossil gas would be tidally driven into 
the primary on a very short timescale near the final merger, 
leading to a bright electromagnetic precursor. Finally in ^ 
we discuss some of the caveats of this model and directions 
for future work. 



2 BASIC PICTURE 

We illustrate the basic scenario schematically in Figure [1] 
A secondary BH with mass Msec and a coplanar accretion 
disc orbits a more massive primary BH at a radius raec- We 
expect the disc to be aligned to the or bital plane of t he two 
BHs because of inclination damping. Ilvanov et al.1 (|l999l ) 
showed that in the case of a BH orbiting a more massive 
BH with an associated accretion disc, the inner accretion 
disc aligns with the orbit of the binary independent of the 
viscosity. As a result, the accretion disc is truncated into two 
parts: an inner accretion disc which surrounds the primary 
and a circumbinary disc which surrounds the binary. Flow 
between the inner disc and outer disc is prevented by the 
presence of the secondary BH at least in one dimension. The 
inner accretion disc has an inner radius, which is typically 
the innermost stable circular orbit (ISCO), and an outer ra- 
dius where it is truncated due to tidal interactions with the 
secondary BH. Similarly, the outer disc is also truncated. 

A similar model has previously been studied by 
lArmitage fc NataraianI (|2002l ). However, we will include sev- 
eral improvements in our treatment. First we include more 
realistic physics for the evolution of the viscous disc, i.e., 
we use a self-consistent alpha disc instead of a fixed viscos- 
ity law. This significantly changes how the disc drains and 
we believe that it provides a more realistic estimate of the 
luminosity evolution of the disc prior to merger. Secondly, 
we focus on the luminosity of the late time electromagnetic 
precursor and its simple shape, which might aid its identifi- 
cation. 



3 BASIC EQUATIONS 

We now present the governing equatio ns for our model. 
These are standa rd in the literature (iLin fc Papaloizmi 
1979bl l 3 , IT986; Ho urigan fc WardI 1 19841: I Ward fc Hourigan 
19891: iRafikovi i2002l : IChand I2OO8I : Ivan de Ven fc Chang 
20081 ) ■ but we st ate them here for cor npleteness. We adopt 
the notation of Icii an^ (|2008h and Ivan de Ven fc Chans! 
(2008,) . We also include the effect of angular momentum 
losses due to GWs, which become important at small bi- 
nary separation. 

The evolution of a viscous disc under the influence of an 
external torque is given by equations for continuity and an- 
gu lar momentum con servation[f| The equation of continuity 
is (|Frank et al.ll2002l ) 

f + (1) 
Ot r or 



where r is the radial coordinate and Vr is the radial com- 
ponent of the v elocity. The equation of ang ular momentum 
conservation is (|Pringlelll98"ll : lRafikovll2002 



a(Sr^n) 1 d{rVrT.r^Q) 
dt r dr 



1 f dT.i. 
' 2-Kr V 



dr 



where 57 = ^/GMBHi?)Jr^ is the orbital frequency, Mbh is 
the mass of the primary BH, and 

^ For simplicity, we do not make a distinction between the center 
of mass of the binary and the location of the primary BH in our 
analysis, even though the secondary is not much less massive than 
the primary in some of our models. 
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Figure 1. Schematic diagram of the model considered in this paper. An accretion disc and a secondary (coplanar) BH with semimajor 
axis rsec surrounds a central BH. The secondary BH opens up a gap in the disc, dividing the disc into an inner disc and an outer disc. 
The inner disc steadily drains into the central BH, while the outer disc viscously spreads and tidally interacts with the secondary. The 
binary's semimajor axis decreases due to tidal interaction with the viscous disc until rgec < ''GWi which point GW emission becomes 
the primary mechanism by which angular momentum is lost. For rscc <S ''GWi the inner disc can no longer viscously respond on the 
merger timescale and is driven into the primary by tidal interactions with the secondary. 



Jvisc = —Inr vl^—-, 
or 



(3) 



is the viscous torque, where v is the viscosity. We have 
adopted the standard a prescription for the viscosity 
jShakura fc Svunvaev|[l97i ) : 



2QPg 



(4) 



with proportionality to Pg, the midplane gas pressure. We 
chose to rely on the so-called "/3 model", i.e., stress propor- 
tional to gas pressure as opposed to total pressure, to ensure 
viscous and thermal stability in our 1-d model. For most of 
our calculations, we use a = 0.10 A simple a-prescription 
applied to radiation pressure d ominated discs is both vis- 
cously and thermally unstabl e jLightman fc EardlevI 1 19741 : 
IPiranI 1 19781 ) . Recent work by iHirose et al.l (|2009l l suggests 
that the thermal instability is avoided in radiation pres- 
sure dominated situations because stress fluctuations lead 
the associated pressure fluctuations. However, the viscous 
instability of radiation pressure dominated discs might re- 
main. Adopting the "/3 model" for the viscous stress will 
only change the evolution at high Mcxt and small radius 
and does not affect our key results: once the binary is at 



3 The exception is the suite of calculations presented in Table 1, 
where we have take a = 0.01 and a = 0.1 



small separation, the discs' evolution is dominated by GW 
losses and not by internal viscosity. 

Using equation ((l|, we simplify ((2} to be 



1 /^Tvis 



dt 



dr 



2nr \ dr 



Solving for Vr in equation ([5} and plugging it into ([T]), we 
find 



9E _ _i_d_ ( djr'^n) 
dt 2-Kr dr \ dr 



(Tvisc — Td) 



(6) 



The orbiting secondary BH exe rts a torque density on 
the disc, dT^/dr, of the fo rm ([Goldreich fc Tremaind 



198d: iLin fc Papaloizoul Il979bl: IWardI Il997l : 



Armitage fc Nataraianll200a : IChandl2008l ') 



or 9 S12 (r - rsec) 



where / is a constant which is 
lArmitage fc NataraianI ([2002). 

Plugging equations Q and ((T)) into we find 



(7) 

0.1 following 



9E 
dt 



ld_ 

r dr 



3r 



1/2 



d_ 

dr 



-/sgn (r - rsec) Q.r'^q^T. I 
97r \{r~rscc), 



(8) 
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where q = Msoc/A/bh is the mass ratio between the sec- 
ondary and primary BHs. 

The radial motion of the secondary BH, assuming a 
circular orbit, is governed by tidal torquing of the disc and 
the torque due to GWs: 



-MsecU 



seci ''sec ' sec 



dt 



— 2gw- 



(9) 



where S — sgn (r — rscc) is the sign of the torque, an d Tpw 
is th e torque exerted by gravitational wave losses l|Petersl 



Taw = -^Mbhc^ 

5 \ Ts, 



7/2 



f Ms, 



Msec + Mbh 



Mb 



,(10) 



where rg = 2GMbh/c^ is the Schwarzschild radius. 

Finally we assume vertical energy balance at every ra- 
dius between heating and radiative cooling, i.e. F = D{r), 
where F is the local radiative flux and D is the local dis- 
sipation rate. In a one-zone model, the local radiative flux 
is 



F ■ 



4 aT^ 
3 kT, ' 



(11) 



where T is the central temperature and k is the opacity, 
which we assume to be constant (Thomson) for simplicity. 
More detailed models would include variations in opacity 
with radius in the disc. We assume that the disc emits as a 
blackbody. The heating is given by a combination of viscous 
dissipation and tidal dissi pation. Viscous dis sipation is given 



dissipation and tidal dissi pation. Viscous dis s 
by the standard formula (|Frank et al.ll2003 ) 



9 ^GAfBH 



(12) 



where Dvisc is the local viscous dissipation rate per unit area. 

We now consider the energy dissipation from tidal 
torques. We assume that the torque on the disc raised by 
the satellite, Td, is mediated by the excitation of spiral 
density waves and locally damped. The timescale associ- 
ated with this tidal torque is the angular momentum of 
the satellite, iscc = -A^sGcf^sccf'sGc, divided by the torque or 
^tido ~ Lscc/Td- The power associated with such a torque is 
then 



Escc 



(13) 



where -Bsec = GMBuMgcc/rsec is the orbital energy of the 
secondary. Because we assume local damping of spiral den- 
sity waves, the corresponding binding energy liberated lo- 
cally (assuming circular orbits) is in proportion to the spiral 
density waves which mediate this interaction. Hence the lo- 
cal tidal dissipation rate is given by 

'dTd/dr\ 



2nr J dr\dTd/dr 



(14) 



As we will show, these tidal interactions are especially 
important at late times when GW losses drive the secondary 
inward on a circular orbit. A tiny fraction {Md^in/Mscc ~ 
10^^", where Md,in = Stt T,r^dr is the mass of the in- 
ner disc) of the angular momentum losses is offset by a gain 
from the tidal interaction with the inner disc. The tidal-GW 
evolution allows this fraction of the prodigious gravitational 



luminosity Lqw < lO^'' ergss~^Qto be released electromag- 
netically by heating the inner disc gas, which is enough to 
power the inner accretion disc close to its Eddington lumi- 
nosity. 



4 RESULTS 

We solve equations (1 8)1 and ([9|l using st andard explicit finite- 
difference methods I Press et al.lll993 ). We choose 200 grid 
points logarithmically spaced between r = risco = 3rg (i.e.. 



the ISCO) and 10 rg. We set a zero-torque boundary condi- 
tion at the inner radius and impose an outer boundary con- 
dition such that there is a constant external feeding rate, 
J\fcxt. As a test, we initially solve these equations for no 
satellite and with inflow boundary conditions and find that 
we recover the steady state a-disc solution with an error of 
< 1%- 

The first two cases we consider are a 10'' M© BH with a 
g = 0.1 and a. q = 0.3 secondary, respectively. We start the 
secondary at an initial radius of rs,o = 10* fg (« 0.01 pc). 
We begin with a low mass (IO'^Mq) disc which extends from 
the ISCO to router = lO^rg. We clear a region around the 
secondary's initial radius, 0.5rs,o < r < 2rs,o, to model the 
initial clearing of a gap around the secondary BH. At the 
outer edge of the grid, we consider an outer Moxt of either 
Moxt = 0.1 Mq yr"^ or 0.01 M© yr'V 

We choose a low initial disc mass (10^ Afe, which corre- 
sponds to a very low accretion rate in steady-state) so that 
our models can be integrated on a reasonable timescale. The 
reason is that the Courant condition, which is set by the vis- 
cous time of the disc near the ISCO, limits the timestep to 
very low values unless the disc has a very low mass. Fortu- 
nately, our results do not depend on this initial condition. 
The mass of the inner disc, Afd,in, approaches an asymp- 
totic value which is independent of its initial value. With 
larger initial disc masses, the inner disc would have initially 
drained on a faster time scale, but the final disc mass would 
remain the same. 

We illustrate this point in Figure[2]where we show mod- 
els with different disc masses of 10^ 10^ and IO^Mq for 
two different mass accretion rates: Moxt = 10~^ MQyr~^ 
(top set of curves) and Afoxt = lO"'^ MQyr~^ (bottom set 
of curves). Note that the mass of the inner disc approaches 
an asymptotic value that is independent of its initial value. 
The reason for this is that the inner disc drains until its vis- 
cous time is comparable to the merger time. At late times, 
however, when the evolution of the system is dominated by 
GWs, the inner disc can no longer viscously drain on the 
timescale of the merger. Hence, the inner disc mass is es- 
sentially frozen. The frozen mass of the inner disc at late 
times is given by the mass of the inner disc when the sys- 
tem transitions from viscosity-dominated evolution to GW- 
dominated evolution, which occurs typically at rpw ~ 500rg 
jArmitage fc Nataraianll2002l : [Haiman et al.ll200^ . also see 
§4.2). In the remainder of this work, we use this fact and 
start with a low disc mass to save computational cost. 



The GW luminosity can maximally be c^/G ' 



ergs s 



but the GW luminosity when the disc is still present is closer to 
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Figure 2. Evolution of the disc mass inside the secondary's or- 
bital radius, rgec, for a 10^ Mq primary and 10^ Mq secondary 
and external mass accretion rates of 10~^ MQyr~^ (Mcxt ~ 
■^Edd; top set of curves), and 10"-^ MQyr~^ (Moxt ~ O.lM^dd't 
bottom set of curves). Three scenarios in each case, with ini- 
tial disc masses of 10-*, 10^, and 10^ Mq, are represented by the 
curves from top to bottom. Note that the asymptotic disc mass 
is independent of the initial disc mass. 



4.1 Viscous Evolution of the BBH and the Outer 
Disk 

We first consider the case of a 10^ Mq BH with a W^Mq sec- 
ondary and an outer mass inflow rate of Mext = O.IM© yr~^ 
(Figure El or Mcxt = O.OIMq yr'^ (Fi gureQ. Figure |3] and 
m show four different cuts from our one- dimensional simu- 
lations to illustrate the different regimes through which the 
binary and the disc evolve with time. In each cut, we plot 
E (solid line matched to the left axis) and fvisc (dashed line 
matched to the right axis). We also show for each cut, the 
merger time of the BBH tmcrgc as a horizontal dotted line 
matched to the right axis. In (a), the satellite and binary 
have just been initialized and the outer disc has filled due 
to mass inflow from the outer boundary. In (b), the sec- 
ondary's radius has decreased by an order of magnitude due 
to its tidal interaction with the viscously evolving outer disc. 
Meanwhile, the inner disc has drained substantially. In (c), 
the binary has shrunk by another order of magnitude in ra- 
dius due to both tidal-viscous evolution and GW losses, with 
the latter becoming dominant at small radii. Note that the 
inner disc density profile has begun to invert, i.e., the sur- 
face density rises with radius as opposed to falling. This is 
due to the inability of the inner disc to viscously respond to 
the increasingly rapid infall of the secondary. By (d) , a large 
surface density spike has built up near the outer edge of the 
inner disc. This forcing of the inner disc onto the primary 
BH causes it to brighten with a maximum luminosity, Lpeak, 
reached at a time, tpeak, before the final merger. 

A crucial aspect of the evolutionary sequence shown in 
Figure [3] and [4] is that, during the early phases (plots a and 



b; before GWs dominate), the viscous time, t^isc = r^/u 
(dashed lines) near the outer edge of the inner disc ad- 
justs to become comparable to the evolving merger time, 
tmcrgc ~ |''s/i'scc| (showu as a horizontal dotted line in plots 
a and b). This implies that the evolutionaxy time for the 
inner disc is the same as the merger time of the binary. This 
adjustment occurs because a more massive disc would drain 
much more rapidly until it drains until the point where its 
viscous time equals the lifetime of the system. Similarly, a 
less massive disc would remain static until the external forc- 
ing from the secondary decreases the disc's outer radius and 
increases the local surface density to the point where the lo- 
cal viscous time matches the merger time. Hence the inner 
disc approaches an asymptotic state that just depends on the 
merger time of the binary. The binary merger time itself de- 
pends on the external feeding rate until gravitational waves 
take over. This self-regulated evolution explains lack of any 
difference of the inner disc mass on the initial disc mass 
shown in Figure [2l Note that once GWs dominate (plots c 
and d), tvisc does not adjust to fmorgc, the disc can no longer 
respond viscously to the rapidly infalling secondary and its 
dynamics becomes tidally dominated as shown vividly in 
panels (d) of Figures [3] and |4l The dissipation of tidal en- 
ergy and the forcing of disc material to smaller and smaller 
radii is the primary means by which energy is produced in 
this late stage. 

As Figure [3] and U show, the general evolutionary se- 
quence of a BBH system can be summarized as follows. Ini- 
tially the BBH starts out at large radii, where the tidal- 
viscous interaction with the outer gas disc shrinks the semi- 
major axis of the BBH. Meanwhile, the inner disc contin- 
ually drains. Eventually, the BBH shrinks until r < raw, 
where GWs become the principle mechanism by which an- 
gular momentum is lost. Shortly afterwards, the inner disc 
can no longer respond viscously to the inspiralling secondary 
and it reaches its asymptotic mass, MasymQ The inner disc 
is forced to smaller and smaller radii by the inspiraling sec- 
ondary until it completely accretes onto the central BH. At 
this time, tpcak before the final merger, the inner disc has a 
peak luminosity, Lpeak- We have run a suite of models at dif- 
ferent a, A/oxt, Mbh, and q to explore the parameter space 
of Md,in, Lpoak audtpoak- Thcsc rcsults are summarized in 
Table 1 and discussed below. 

4.2 The Electromagnetic Precursor of the BBH 
Merger 

As Table 1 indicates, the peak luminosity of the electromag- 
netic precursor can be appreciable even though the asymp- 
totic mass of the inner disc seems insignificant. The pre- 
cursor's light curve is shown as a function of time prior to 
merger in Figure [S] where we plot the luminosity due to 
tidal and viscous dissipation as a function of the time be- 
fore merger, tminua (measured from when the BHs finally 
merge, which is not the same as the instantaneous merger 
time raec/|usec|) for the two cases of Mbh = lO^Af© (a) and 
10* M0 (b). The results shown in Figure [S^ are for the two 
mass infiow rates, O.IM© yr~^ (solid line) and O.OIM0 yr~^ 

^ For consistency, we take the asymptotic disc mass, Maaym, to 
be Af(j i„ at r = lOOrg. 
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Figure 3. Snapshots of the disc properties for a IO^A/q primary BH with a 10® Mq secondary BH (whose position is marked with 
an "x"), for an external accretion rate of 10~^MQyr~^. We plot the surface density E (solid lines), viscous time tvisc (dashed line, 
matched to right axis), and merger time of the secondary tmorge (dotted horizontal line, match to the right hand side), (a) shortly 
after initialization, when the outer disc has accumulated some mass due to accretion, (b) after the semimajor axis has decreased by an 
order of magnitude due to tidal-viscous interaction with the outer circumbinary disc. Note that the mass of the inner disc has declined 
significantly compared to (a) while the outer disc is approximately an alpha disc with a cutoff around the radial position of the secondary, 
(c) the semimajor axis has decreased by yet another order of magnitude, and the inner disc mass has reached it asymptotic value. Angular 
momentum losses are now dominated by GWs. (d) the secondary tidally forces the remaining mass in the inner disc to accrete rapidly 
into the central BH. It is this final forced feeding of the central BH which leads to the electromagnetic precursor (see Figure[5K,b). 



(dashed line), and for q = 0.1 and q = 0.3 (top and bottom 
plots respectively) . The electromagnetic precursor manifests 
itself as a sudden brightening to « 0.1 — 1 1/Edd during the 
last few days before merger. There are a few notable char- 
acteristics of the hght curve. The peak luminosity depends 
only on the mass of the inner disc. The time of the peak 



ipcak depends only on the masses of the two BHs. Finally, 
the light curve of the precursor has a very characteristic 
power-law shape. 

The first two points are clear from examining Figure 
\5ii,h. For instance for the Mbh = lO'^M© {10^ Mq) case in 
Figure [5^(b), the characteristic shape and time of the peak 
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Figure 4. The same as Figurc|3] but for an external accretion rate of 10 ^Mq yr ^ . The main notable difference is that the asymptotic 
mass of the inner disc is significantly smaller than for Afcxt = 10~^ Mq yr~^ case of Figure[3] 



of the light curve are the same for the two different Moxts, 
but the normaUzation is different. The normalization of the 
light curve depends on the mass of the inner disc, which in 
turns depends on the outer mass accretion rate and a (see 
Table 1). On the other hand, the timescale of maximum 
luminosity only depends on the masses of the two BHs as 
one can see by comparing the cases g = 0.1 and q = 0.3 in 
Figure [5^,b. 

The light curve follows a characteristic power law shape, 
which we now derive. Tidal interactions between the sec- 
ondary BH and the accretion disc drive the outer edge of 
the disc to smaller radii and the dissipation of the induced 



density waves leads to heating of the disc. The amount of 
energy dissipated at any time is 



GManiWd 1 

''merge ; 



(15) 



where is fixed as we are in a regime with the merging 
time much shorter than the viscous time. This mass is set 
primarily by the viscosity of the disc, i.e., a, and the external 
Mext . It directly influences the luminosity of the inner disc in 
this GW driven phase of the merger. The merging timescale 
due to GWs is, 
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Table 1. Asymptotic inner disc mass and several quantities charac- 
terizing the EM precursor for a range of black hole masses, accretion 



rates, and a. We show Mbh 
[Mq], tpeak [days], and Lp;,^^ 


[Mq], q, 
[ergs s ■'■] 


Moxt [Mq yr-1], a 


Mbh q A/oxt 


a 




^pcak 


-^pcak 


10 




X iu 




X 10 


0.1 0.1 


0.1 


1.2 


0.07 


1.9 




0.01 


12 


0.07 


3.6 


0.01 


0.1 


0.12 


0.07 


0.2 




0.01 


1.2 


0.07 


0.38 


0.3 0.1 


0.1 


0.1 


0.04 


0.5 




0.01 


1 


0.04 


0.3 


0.01 


0.1 


0.008 


0.04 


0.05 




0.01 


0.1 


0.04 


0.027 


1 n7 
10 




X iU 




X 10 


0.1 0.1 


0.1 


6.6 


0.7 


2.0 




0.01 


460 


0.7 


14 


0.01 


0.1 


0.73 


0.7 


0.23 




0.01 


6.4 


0.7 


2.0 


0.3 0.1 


0.1 


0.62 


0.4 


0.19 




0.01 


5.6 


0.4 


1.5 


0.01 


0.1 


0.062 


0.4 


0.019 




0.01 


0.75 


0.4 


0.2 


10« 




xlO-2 




xl0« 


0.1 1 


0.1 


20 


7 


6.5 




0.01 


90 


7 


28 


0.1 


0.1 


2.9 


7 


0.9 




0.01 


18 


7 


5.7 


0.3 1 


0.1 


2.4 


4 


0.67 




0.01 


15 


4 


4.2 


0.1 


0.1 


0.33 


4 


0.09 




0.01 


3.3 


4 


0.9 


^ 5 a Mbh { a 

Imcrgc , , 1 

8 C Msec V '^9 








( 



(16) 



Combining equations (|15p and (|16p . we find 

L cc t^lt- (17) 

This power law shape is a good approximation to the hght 
curve in the more detailed calculations in Figure [5K,b. 

The associated effective temperature is T^h oc 
{L/rd{t)^y^'^ oc t^Jl^^ ■ Figure [U] shows the corresponding 
multi-temperature blackbody spectra at tminus = 1 day and 
10 days for Mbh = lO^M© (left panel) and 10* M© (right 
panel). These spectra are also shown for different mass ra- 
tios, q = 0.1 (thick lines) and q = 0.3 (thin lines), and for 
different outer accretion rates, Mcxt = MEdd (solid lines) 
and O.lMsdd (dashed lines). As this plot demonstrates, the 
majority of the emission prior to merger will be in the ex- 
treme UV and soft X-ray. As a result, the electromagnetic 
precursor may be difficult to detect from the ground or in 
the presence of obscuration. We have presumed no reprocess- 
ing of the inner disc emission by the outer disc. If present, 
this may increase the luminosity of the source in more easily 
detectable wavelength, i.e., optical or near-IR, with possibly 
significant time delays. Large amplitude spiral density waves 
from the strong tidal forcing may contribute to non-thermal 
high energy emission. One notable aspect of these spectra 



is a gap in the emission between the low energy part, which 
is due to emission from the circumbinary disc, and the high 
energy part, which is due to emission from the tidally forced 
inner disc. 

The mass in the inner disc can be dramatically affected 
if additional sources of mass inflow exist beyond the viscous 
condition considered here. For instance, if the gap between 
the inner disc and outer disc is not completely clear, then 
a small trickle of mass from the outer disc to the inner disc 
could change the asymptotic inner disc mass significantly. To 
quantify this point, at the transition radius between viscous 
and GW driven evolution, the mass drainage rate of the 
inner disc is 

. ^10-^°-10-^Meyr-\ (18) 
imcrgc [raw ) 

where fmorgc = ^scc/wsoc is the instantaneous merger time. 
This very small mass inflow rate implies that any addi- 
tional source of mass whose time averaged input rate over 
tmorgc ~ 10®~^yrs equals or exceeds equation (|18|l would 
greatly modify the inner disc evolution. The calculated 
asymptotic disc mass may thus be considered to be a con- 
servative lower limit on the mass of the inner disc. 



5 DISCUSSION AND CONCLUSIONS 

We have calculated models of BBH evolution subject to 
torques from a circumbinary accretion disc and gravitational 
radiation losses. In general, the viscous time inside the orbit 
of the secondary BH is much less than the timescale for tidal- 
viscous interaction with the circumbinary disc to shrink the 
orbit of the binary. At flrst glance, this suggests that the 
inner disc should completely drain away. However, we have 
shown that the inner accretion disc maintains a small, but 
nontrivial residual mass (Fig[2]), the value of which is set 
when GWs take over (at row ~ 500rg) as the primary 
mechanism shrinking the orbit of the binary. This inner 
disc, which is then tidally forced into the primary by the 
secondary as it spirals inward, brightens significantly dur- 
ing the last day prior to merger. The bolometric luminosity 
rises as t~fn*a, peaking at one day with a peak luminosity of 
~ O.lLsdd for a representative lO'^ M© BH (see Table 1 for 
results for various as, Mbh, Mcxt, and q) 

Such an electromagnetic signature is very important for 
the identification of the host galaxies of BBH mergers, but 
our study is subject to several caveats. One such caveat is 
related to the persistence of a thin accretion disc at very low 
accretion rates. Studies of X-ray binaries suggest that below 
an accretion rate of lO"** — lO'^A^Edd, accretion discs can 
transition from a thin radiatively efficient accretion disc to a 
thick radiatively-in efficient accretion flow (iGallo et al.ll2003 : 
i Fender et al.ll2004l ) . If this were the case in our scenario prior 
to GW driven evolution, the inner disc would completely 
drain into the primary BH and hence the characteristic elec- 
tromagnetic signature we have identified would be absent. 
We emphasize, however, that the inner disc asymptotically 
approaches a gas-supported, optically-thick thin disc solu- 
tion, which in itself is stable at very low accretion rates. 
Thus it is not apriori obvious whether the transition to a 
geometrically thick disc must occur in all cases. 

Our calculated mass for the inner disc represents a con- 
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Figure 5. Bolometric luminosity (thick curves) and mass of the inner disc (thin curves) as a function of time before merger, tminua- 
(a) Evolution for a IO^Mq primary with a 10^ Mq secondary (upper panel) and a 3 X 10^ Mq secondary (lower panel). We plot 
this evolution for Mext = lO~^M0yr~^ (solid lines) and Mext = lO~'^M0yr~'^ (dashed lines), (b) Similar evolution for a 10* 
primary with a lO^M© secondary (upper panel) and a 3 X IO^Mq secondary (lower panel) and for Mext = 1 Mq yr~^ (solid lines) and 
A-/cxt = 10~^ Mq yr"-*^ (dashed lines). Note that the precursor approximately follows a t„,j'j{ug power law until the peak. 




Figure 6. Multi-temperature blackbody spectra at tminus = 1 da-y for Mbh = Mq (left) and at tminus = 10 days for Mbh = IO^Mq 
(right) in the evolutionary scenarios of Figure [5] Thick lines represent a mass ratio q = 0.1 and thin lines represent q = 0.3. Note the 
gap in the emission between the low energy component, which arises from the circumbinary disc, and the high energy component, which 
arises from the tidally forced inner disc. 
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servative lower limit. When the secondary transitions from 
a viscosity driven evolution to a gravitational wave driven 
evolution, the inner disc mass corresponds to a mass accre- 
tion rate of M ~ md/tmerge ~ 10~® Mq yr~^. For any time 
averaged mass source of cold gas present over tmagc ~ 10®- 
10^ yrs that is greater than this, the mass of the inner disc 
would be larger than what we have calculated. For instance 
the mass loss from a few massive s tars such as found in 
the Galactic center ijOuataertl l2004l : iLoeh ^20 04) or an old 
popu lation of stars as in the nucleus of M31 (|Chang et al] 
I2OO7I ). could easily enhance the mass of the inner disc. In 
addition, our one dimensional models presume that the gap 
between the inner and circumbinary disc is clean and that no 
mass flows between the two. While this may not be an un- 
reasonable assumption given the comparable masses of the 
two BHs, two-dimensional simulations of gaps induced by 
protoplanets in a protoplanetary disc suggest that ~ 10% of 
Mext in the outer disc can flow across the gap onto the inner 
disc (iLubow fc D'Angeioll2006l : [MacFadven fc Milosavlievid 
l2008l 'l. Mass leakage across the gap in our case is likely to 
be much smaller than what it is in the protoplanetary case, 
but even a small amount of mass will make a substantial 
difference to the evolution of the inner disc. 

We expect the perturbation to the GW signal from the 
tidal forcing of the inner disc to be minor in the sense that 
it would not be detectable with LISA. In the final year 
of inspiral, the inner disc has reached its asymptotic mass 
of ~ 10^^ — 1O~*M0 (see Figure 5a,b). In this stage, the 
inner disc modifies the torque from GWs by a factor of 
Md/Msec ~ 10"^'^. Such a tiny mass correction is well below 
the best achievable precision on mass meas urements for the 
secon dary with LISA (e.g., iHughea ,200 j : iLang fc Hughes! 
l2008l ). 

The tidal forcing of an accretion disc by a spiraling bi- 
nary is a somewhat unusual situation in astrophysics and 
the possibility that it occurs in the mergers of BBHs mo- 
tivates further exploration of the physics of this scenario. 
One interesting possibility is that the precursor re bright- 
ening event that we have identified could be time-variable 
due to perturbations from the orbiting secondary. This is 
analogous to the case of superhumps, where perturbations 
from a low mass secondary drives the disc in a cat aclysmic 
variable system to be mildly eccentric (|Lubowlll99lh . A sim- 
ilar eccentricity growth ha s been su ggested for the outer 
disc of BBHs by M acFa dven fc Mil osavlievid l|2008l ') and 
ICuadra et al.l l|2009l ). Time variability on the binary's pe- 
riod (or some multiple) would allow for an unambiguous 
association of the electromagnetic signal with the GW sig- 
nal measured from future gravitational wave detectors such 
as LISA. The emergence of a time variable signal from the 
tidally forced accretion of this fossil gas is worthy of further 
investigation. 
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